Ab initio investigation of VOSeC>3, a spin gap system with coupled spin dimers 
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Motivated by an early experimental study of VOSeOs, 
which suggested that it is a quasi-2D system of weakly cou- 
pled spin dimers with a small spin gap, we have investigated 
the electronic structure of this material via density-functional 
calculations. These ab initio results indicate that the system 
is better thought of as an alternating spin- 1/2 chain with 
moderate interchain interactions, an analog of (VO)2P207. 
The potential interest of this system for studies in high mag- 
netic field given the presumably small value of the spin gap 
is emphasized. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



The search for low-dimensional s=i materials with a 
singlet-triplet gap in the spin excitation spectrum has 
been a very active field of research since the discovery of 
high T c superconducting cuprates. The original motiva- 
tion came from the possible connection between spin-gap 
and superconductivity, an issue still intensely debated in 
the context of cuprates. It became clear however with the 
extensive investigation of spin ladders that even purely 
insulating magnets with a spin gap can exhibit a very 
rich physicsEra. The basic idea is quite simpleO: When 
the gap of a system consisting of weakly coupled dimers 
is closed by a strong enough magnetic field, some dimers 
are promoted to triplets. These triplets behave as a 
gas of quantum particles, and they can undergo differ- 
ent kinds of phase transitions as a function of magnetic 
field or temperature depending on dimensionality, topol- 
ogy, etc... To reach this regime however, one needs sys- 
tems with moderate gaps - up to 20 K or so. This is 
the typical range of magnetic coupling in organic mate- 
rials, and the first investigations have been accordingly 
carried out on organic systems. These systems are tricky 
however, and even the most extensively studied organic 
ladder, C^^sHwij^^CU, is still a subject of debate 
and controversiesHa. The difficulties with organic sys- 
tems are two-fold: the exchange paths are not directly 
evident on the basis of the structure, and it is in most 
cases impossible to grow large single crystals. 

These problems are usually much easier to solve in 
inorganic transition-metal compounds. The main diffi- 
culty there comes from the order of magnitude of the 
exchange integrals - typically several hundreds of degrees 
Kelvin - but if the geometry is such that the exchange 
integrals are reduced substantially in the spirit of the 
Goodenough-Kanamori-Anderson rulesQ, then the low- 
temperature properties are usually easier to reach. Sev- 



eral qualitatively different classes of behaviour have been 
identified so far in this family, among which are incom- 
mensurate spin-spin correlations in the spin-Peierls high- 
field phase of CuGeOala, magnetic field- induced long- 
range order in TICUCI3 and KCuCl30i3, and magneti- 
zation plateaux in SrCu2(B03)2tJ due to magnetic field- 
induced localization of spin-triplets. The behaviours re- 
ported so far do not exhaust the possibilities though. 

In that respect, the vanadyl system VOSe03 is poten- 
tially very interesting. This compound was first investi- 
gated by Trombe et aZE3, who synthetized the system, 
determined its structure, and measured the temperature 
dependence of the susceptibility. They concluded that 
the system consists of weakly coupled dimers, and that 
the residual couplings are presumably of two-dimensional 
character. Besides, although the susceptibility data re- 
ported by Trombe et al. do not allow the determination 
of the spin gap because of a large Curie tail, the location 
of the maximum of the susceptibility suggests that the 
gap is probably in the appropriate range to be closed by 
a magnetic field. However, the details of the dominant 
interactions, their origin as well as the role of the active 
orbitals at the Fermi surface could not be quantitatively 
discussed in that work, and the conclusions are only pre- 
liminary. 

In order to check the validity of this analysis, and in 
particular to determine whether the couplings are indeed 
essentially two-dimensional, we have performed an ab ini- 
tio calculation of the electronic properties of this system, 
followed by a tight-binding-downfolding analysis in order 
to define the important hopping parameters in this com- 
pound. As we shall see, the emerging picture is somewhat 
different, although. not less interesting, than that guessed 
by Trombe et alS3. 

In order to understand the electronic behavior of this 
system, an examination of the crystal structure is es- 
sential. VOSeOs crystallizes in the monoclinic space- 
group P2i/c with lattice parameters a = 4.0168 A, b = 
9.788 A, c= 8.001 A, (3 = 99.42°, and it contains four 
formula units per primitive unit cel£3. The vanadium 
ions V + form square pyramids with the neighboring oxy- 
gens. Two adjacent up and down square pyramids share 
an edge and build dimeric [V^Os] 8- units. These units 
form chains along the x direction (see Fig. [j] (a)) and 
are linked through Se atoms in the yz plane as shown 
in Fig. [I] (b). The vanadium atom in the square pyra- 
mid is shifted towards the apex oxygen and forms a short 
vanadium-oxygen bond (d = 1.611A) characteristic of a 
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vanadyl ion V0 2+ . The distance between two V 4+ ions 
within the dimeric unit is of d = 3.175A while the dis- 
tance between two dimeric units along x is d = 4.015A 
and along z is d = 4.841A 

We have carried out a first-principles study based on 
the density-functional theory (DFT) in order to derive 
the electronic properties of VOSeC>3. We haise used 
the generalized gradient approximation (GGA)Ej in or- 
der to include the non-local effects within the gradient 
approximation to go beyond the local density approxi- 
mation (LDA). Calculations have been performed within 
the framework of both the full-potential linearized aufc 
mented plane wave (LAPW) method based on W1EN97 — 
code and the linearized muffin tin orbital (LMTO)t2 
method based on the Stuttgart TBLMTO-47 code. The 
results obtained by both methods are in agreement with 
each other. 

In Fig. H we show a plot of the band structure along 
the symmetry path in the Brillouin zone T — (0, 0, 0), 
B=(-7r,0,0), D=(-7r,0,7r), Z=(0, 0, tt), T, Y=(0,7r,0), A=(- 
tt, tt, 0), E=(-7r, tt, tt). This system shows four narrow 
bands (four since the unit cell has four vanadium atoms) 
at the Fermi level with a bandwith of about 0.5 eV. These 
bands are of vanadium 3d xy character (in the local frame 
of reference with the z axis pointing along the axis con- 
necting the vanadium ion and the apical oxygen of the 
square pyramid) with small mixing of oxygen p charac- 
ters. They are separated by a gap of about 2 eV from the 
lower valence bands and a gap of 0.1 eV from the higher- 
unoccupied conduction bands. The system is half-filled 
and the insulator behavior observed in this system should 
be explained by the effect of electron correlation which is 
not taken fully into account in the LDA or GGA calcu- 
lations. Note that the dispersion along the x axis (chain 
direction) is small while in the yz plane the dispersion is 
of about 200 meV. This behavior already indicates that 
the important interactions in this system will be in the yz 
plane. We shall analyze this point in more detail in terms 
of an effective hamiltonian with parameterized hopping 
integrals. _ 

We have employed LMTO-based downfoldingEi and 
tight-binding analysis on the ah initio results to obtain 
an effective few-orbital tight-binding description of this 
material. The downfolding method consists in deriving 
a few-orbital effective Hamiltonian from the full LDA or 
GGA Hamiltonian by downfolding the inactive orbitals 
in the tails of the active orbitals kept in the basis chosen 
to describe the low-energy physics. This process results 
in renormalized effective interactions between the active 
orbitals, V d xy in the present case. By Fourier trans- 
forming the downfolded hamiltonian — ► Hr, 

H R = -Y d t ij{c]c i + c\c j ) (1) 

one can then extract the effective hopping matrix ele- 
ments tij between the vanadium ions. The comparison 
between downfolded-tight-binding bands and the DFT 



bands are shown in Fig. g. The band dispersion of the 
four-band complex close to the Fermi level can be re- 
produced well by considering a few short-ranged hop- 
ping parameters (see Fig. td is the intradimer hop- 
ping between two edge-sharing V 4+ Os pyramids and it 
is expected to give the largest contribution since itpCprre- 
sponds to a superexchange path V-O-V (a ~ I03°)Ej. Its 
value is essentially related to the width of the four-band 
set. In Fig. |^ we show the projection of the electron den- 
sity on the zy plane (compare with Fig. |f](b) ) , where the 
contribution of the Vd^-Op-Vdzy path to the intradimer 
coupling can be clearly identified. 

The next important hopping integral is t<i which des- 
cribes the interdimer coupling along the z direction. 
Other hopping parameters in the yz plane, ti, t\, tj, are 
smaller but non-negligible. They are responsible for the 
band-splitting along TD and dispersion along ZY. Note 
that the pathes described by t^ and t± are not equivalent 
as can be observed in the electron density plot of Fig. |] 
(compare with Fig. The hopping integral along the 
chain direction x, t x proves to be very small, which in- 
dicates that this interaction path is negligible. This is 
to be expected from the fact that the active d^j, orbital 
defined in the rotated local frame of reference with the z 
axis turned along the vanadium -apical oxygen, has very 
weak interaction with the apical oxygen p orbital. Still 
there are two more parameters t v and t s which are not to 
be neglected since they are responsible for the band split- 
ting along the path DZ and dispersion along X-A- This 
system shows a-uery similar behavior to CSV2O5E3 as well 
as (VO)2P207lLZl though in those cases important pathes 
of interaction where provided through V 5+ 04 and P 5+ 04 
tctrahedra groups respectively, which are not present in 
VOSeC>3. Instead, Se 4+ C>3 trigonal pyramid groups con- 
tribute decisively to the interdimer interaction pathes in 
VOSeC>3. Those groups can be considered as playing 
the equivalent role to the tetrahedra groups in the above 
mentioned compounds. 

A detailed comparison between the hopping parame- 
ters in VOSeOs and CSV2O5II3 confirms the similitudp_pf 
their behavior (compare td, t%, £4 with t\, t^, t$ in Ref.113). 
There is a quantitative distinction though between both 
systems which is related to the band-splitting in VOSeOs 
along the path DZ described by the hopping parameters 
t v and t s which is not present in CSV2O5. 

The analysis of the effective tight-binding model for 
VOSeOs leads us to conclude that this system shows a 
coupled spin dimer behaviour with important interdimer 
coupling along the z direction and moderate to negli- 
gible couplings along y and x. This picture provides a 
quantitative description of the important inter.actions in 
VOSeOs as opposed to that of Trombe et aZ£3 on the 
basis of their susceptibility data. 

An estimate of the exchange integral related to the 
dominating interdimer interaction parameter td can be 
obtained by using the relation J ~ M^/U where U is the 
effective onsite Coulomb repulsion on the vanadium site. 
Such an estimate is valid since the path described by td 
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corresponds to a V-O-V superexchange path. Values of-U 
~ 4-5 eV have been proposed for other vanadyl systemsliij. 
Assuming that this value is similar for V0Se03 we get a J 
~ 5 meV ~ 55 K which is of the same order of magnitude 
as the J value obtained by Trombe et al. (J ~ 30K) from 
the analysis of their susceptibility data. 

We hope that the present results will motivate further 
experimental efforts to understand the properties of this 
system. It would be particulartly interesting to inves- 
tigate its properties under high magnetic field since its 
topology is a priori different from that of other materials 
in which the spin gap could be closed. 

We acknowledge useful discussions with P. Millet and 
the support of the Swiss National Fund and the Deutsche 
Forschungsgemeinschaft . 
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FIG. 1. Crystal structure of VOSe03 projected in the (a) 
zx plane and (b) yz plane. The large black circles are V 4+ in a 
square pyramid environment of oxygens (small white circles). 
Se atoms are represented by grey balls. Shown in this figure 
are also the relevant hoppings in this material. 
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FIG. 4. Projection on the zy plane of the electron density 
for bands close to the Fermi level (compare with Fig. [l]). Note 
the Vdxy-Op-Vdxy pathes contributing to the intradimer cou- 
pling in VOSe0 3 . 



FIG. 2. Band-structure for VOSe3 along the path 
F-B-D-Z-r-Y-E. The boxes indicate the band character in the 
local coordinate system. 




FIG. 3. Comparison of the rfouin/oWerf-tight-binding bands 
(solid lines) with the DFT bands (dotted lines). The 
tight-binding parameters (see Fig. [lj) are (in eV) t^= 0.083, 
t 2 = 0.079, h= 0.012, t 3 = 0.005, t 4 = 0.040, t v = 0.040, 
t s =0.016, e =-0.034, t x = 0.0. 
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